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ABSTRACT 

In  this  note  we  determine  the  natural  parameter  space  of  the 
iterated  FGM  distribution,  and  thereby  show  that  the  maximum  correlation 
is  higher  than  what  was  previously  known.  We  also  show  that  one  single 
iteration  can  result  in  tripling  the  covariance  for  certain  marginals. 
Finally,  we  show  that  there  exist  no  marginals  for  which  the  single 
iteration  will  bring  about  higher  negative  correlation. 


1. 


INTRODUCTION 


The  bivariate  Farlie-Gumbel-Morgenstem  (FGM)  system  of 
distributions  [see  Johnson  and  Kotz  (1975)]  has  joint  cumulative  distri¬ 
bution  function  of  the  form 

H(x,y)  *  F(x)G(y)[l  -  aF(x) G(y)  ]  ,  (1) 


where  F  and  G  are  the  marginal  cumulatives,  T=l-F,  G=l-G  and 
a  is  a  real  number  so  that  H  is  a  bona  fide  distribution  function. 

The  FGM  distribution  has  a  considerable  appeal  in  model  building 
since  it  provides  a  convenient  way  to  construct  a  joint  distribution 
having  specified  marginals.  For  details  on  applications,  see  Conway 
(19S5)  and  Schucany,  Parr  and  Boyer  (1978).  The  usefulness  of  FGM, 
however,  was  marred  by  the  fact  that  it  is  restricted  to  describe 
relatively  weak  dependence.  As  was  pointed  out  by  Schucany  et  al. 
(1973),  if  X  and  Y  are  bivariate  H  with  absolutely  continuous 
marginals  then  the  correlation  coefficient  p  between  X  and  Y  can 
never  exceed  1/3. 

Various  attempts  to  enlarge  the  family  have  been  made.  Johnson 
and  Kotz  (1977)  proposed  the  following  "iterated"  generalization  of 
(1): 


FG{1  ♦  OjFff  ♦  a2FG?ff  ♦  ...  ♦  ^(FG) 


[4i] 


Clearly,  H  has  marginals  F  and  G  .  We  shall  call  (2)  the  "k-l 
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(7) 


2  "2  7  "7 

Cov(X.Y) -BjCj)  (l'i2-lJ22)(v12"v22) +°2(2}  (u23"W33)  (v23'V33)  +  *  *  *  * 


a1^22"Wll^V22"vll^  +  a2^35‘W22^V33‘V22^  +  ***  » 


(8) 


last  equality  following  from  recursive  formula  for  ejected  values  of 
order  statistics.  Note  that  Cov(X.Y)  depends  on  F  and  G  only  thru 

wnn  and  vnn'  n=1»2» • • • • 

EXAMPLE  1,  Specializing  (3)  to  F=G=  Uniform  (0,1),  we  see  that 


Gov(X,Y)  «  o/36  +  6/144 


(9) 


p  *=  a/3  +  6/ 12 


(10) 


EXAMPLE  2.  Again  for  FGM  with  one  iteration  (3),  let  F  *  G  = 
Normal  (0,1),  un  =  «  0,  P22  “  V22  =  7r"55’  w33“  v33*  ) » 


Cov(X,Y)  *  p  *  it  *a  +(4jt)  *6 


(ID 


2.  NATURAL  PARA^TER  SPACE 

The  natural  parameter  space  0  of  (2)  consists  of  all  (Qi’a2>***» 
ak)  such  that  (2)  represents  a  bona  fide  cumulative.  For  k»l  case,  Q 
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has  been  determined  by  Canbanis  (1977).  In  particular,  for  absolutely 
continuous  F  and  G,  R  =  {a:  |a|  s  1}  . 

For  ks  2  the  structure  of  0  is  conplicated  and  is  difficult 
to  describe  in  closed  forms.  Unfortunately,  the  usefulness  of  the  FGM 
family  (2)  depends  on  how  rich  the  family  is.  In  other  words,  how  big 
is  11? 

In  this  note  we  determine  (2  for  k  =  2  and  for  absolutely 
continuous  F  and  G  . 

1.  For  absolutely  continuous  F  and  G,  the  natural 
parameter  space  R  of  (2)  is  convex.  - 

— :£?„?.•  Let  0<p<l,  =  (a^, —  ,s»j.)  c  H  and  *  (<Xj, . . .  )e  ft  . 

By  definition,  hJx.yJaO,  h^x.yJiO  a.e.  It  folloxvs  easily  that 

kpa+(l-p)u*(x>y)  =  Pha(x*y)  +  Cl-p)hQ*(x,y)  0  -a.e.  D 

For  F  =  G  =  Uniform  (0,1),  let  us  write  the  density  of  (3)  as 

h(x,y)  =  1  ♦  oc(l-2x)(l-2y)  +  Bxy(2-3x)  (2-3y)  ,  0s*sl  ,  (12) 

and  define 

J  2 

A  ■  ( (a , B)  :  M*l,  B sir"!!9'*?'38  .}  .  (13) 
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It  is  easy  to  verify  that 


In  order  to  prove  A  c  (I  let  us  first  observe  that 

(i)  A  is  the  convex  hull  of  A  U  {(1,-2), (-1,0)},  where  A* 
is  the  north-east  boundary  of  A  : 

A*  »  {(a,8)  :4iaS  l  ,  6  »  £*l£**Z jgl }  -  (i4) 

and  (ii) 

(a, 8)  e  A*  a2  -  8(3  -  o)  ♦  $2  -  0  .  (IS) 


In  view  of  the  convexity  of  ft  ,  it  suffices  to  show 


((1,-2), (-1,0)}  c  ft 


(16) 


and 


A  e  ft  .  (17) 


(16)  is  easily  verified.  To  prove  (17)  let  us  note  that  if  (a, 8)  e  A. 
then  h  a  is  nonnegative  on  the  boundary  of  unit  square: 

CtP 


haS(x,y)2°  x-  0  or  x*l  or  y»0  or  y*  1  . 


(18) 


Let  A* 


ys 


u(a)  =  < 


a+S+/c  -»8+32 
33 


0  2  0 


28-a+»42+2o3+432  os0 
68 


(19) 

(20) 


The  y-section  of  h  ,  t  quadratic  function  of  x  ,  is  either  concave  or 
monotone  on  (0,1),  thus  attains  its  minimum  at  the  end  point  x=0  or 
1  ,  which  is  known  to  be  nonegative  by  (18).  Hence,  in  order  to 
ensure  h(x,y)  2  0  we  hereafter  restrict  our  attention  to  y>u(a) 
and,  by  symmetry,  to  x  >  u(a)  . 

For  y>u(a),  min  h(x,y)  is  attained  at 


6 


X  -  X  =  8y(3y-2)-a(2y-l) 

y  "  38y(3y-2) 

2 

Notice  that  for  a*0,  x  s  <  u(a) ,  a  case  already  disposed  of.  Here- 

y  ^ 

after  we  assume  a<0  (thus  u(a)  is  given  by  (20)).  We  leave  the 
reader  to  check,  making  use  of  (15),  the  following: 

-3Sy(3y-2)h(xy,y)  -  9B2y4  -  6B(2B-a)y3 

+  (4a2+432-7a8-93)y2  +  (66+2uB-4a2)y + u2 

*  y(y-l)[(3y-2)Qa(y)  *  a2]  ,  (21) 

where 

Qa(y)  =  9g2y2  *  53(2a+S)y  t  a2  .  .  (22) 

Q  (y)  is  strictly  increasing,  stTictly  positive  on  y>0,  thus 

Ci 

2 

(3y-2)Qa(y)  >  0  and  is  strictly  increasing  on  p  y  .  In  order  to 
show  LHS(21)  <  0  for  y  >  u(a)  it  suffices  to  show 


(3u-2)Qa(u)  -  o2  >  0  .  (u=  u(a)) 


(23) 


But  (20)  =*> 


68u  -  2(23-u)u-  a 


(24) 
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t 


(24)  ♦  (20) 


LHS(23)  =  -  j  a(36u-g+a) 


j  a(a+B+*42+2oB+4B2)  >  0  , 


completing  the  proof  that  Acfi.  .*.  .  -for  -fU  - 

To  prove  A=£I  for  arbitrary  absolutely  continuous  F  and  G  , 

* 

let  us  notice  that  hQe(x,y)  =  f(x)g(y)hog(x,y) ,  where 

* 

h*3(x,y)  *  l+a[l-2F(x) ] [l-2G(y) J  +  8F(x)G(y) [2-3F(x) ][2-3G(y) ]  . 

Thus  it  is  equivalent  to  the  problem  of  showing 

05nSI  • 

reducing  the  problem  to  the  uniform  marginal  case. 

3.  POSITIVE  DEPENDENCE 

What  is  the  maximum  dependence  attainable  in  the  FGM  family  (2)? 
Clearly  if  we  allow  F  and  G  to  be  discrete,  even  in  the  noniterated 
case  (1),  it  is  possible  to  attain  max  p  =  1  . 

For  absolutely  continuous  case,  the  problem  was  studied  by  Kotz 
and  Johnson  (1977)  and  Schucany  et.al.  (1978),  both  assuming  jojJ  s  1  , 
|a2l  s,l,....|ak|  si  and  k  •  . 


We  shall  restrict  our  consideration  to  one  iteration,  absolutely 

continuous  case  (12),  with  finite  a  and  a  , 

X  / 


{hfi0  :  (a,g)en) 


As  was  pointed  out  in  (7),  for  given  F  and  G,  Cov(X,Y)  depends  only 
on  and  ,  which  in  turn  dictates  the  optimal  choice  of  a 

and  0  . 


Cov  (X,Y)  .  o  v22'“l2  u22'v12  S  v33  u23 

”,  °y  '  4  °y  °x  °y 


For  ex-rpies  1  and  2,  the  optimal  choice  of  (a, 6)  is  a  point  on  the 

*  2  2 
bound  ay  ;et  A  with  slope  n)  /(u---^)  »  namely, 

at  =  —  -  -  3,  3  =  2-  -  ,  yield: ru- 


m- 1 


/13-1 


=  .45426  (Example  1) 

=  .41469  (Example  2) 


For  general  case,  to  obtain  a  bound  for  (25),  we  note  that  for  any 


F  with  finite  , 


u22-v 


M33_W23 


(Ludwig  (1973)).  Substituting  (28)  into  (25),  we  get 


P  <  1  a  +  _31  e 

3  240  (29) 


* 

RHS  above  is  maximized  over  A  at 


a  «  129  _  1 

T/4881 


e  =  2 


18  # 
-^4881 


Thus 


Pi 


_l(  1627 
40  l  -^4881 


.5072 


(30) 


In  both  Examples  1  and  2  we  see  that  one  single  iteration  has 
improved  p  from  1/3  to  (/T3-l)/6,  and  from  a”*  to  (/U-l)/(2ir) , 
respectively.  In  both  cases,  the  improvement  is  (/l3-3)/3,  or,  a 
little  over  30%. 

Our  next  example  shows  that  one  single  iteration  could  result  in 
nearly  200%  increase  in  covariance. 


Example  3.  Let  F=G  be  Pareto  truncated  at  t(>l)  : 


f(x)  =  try  (i-  7)  lsxst  . 

2  1 

Then  p1]L  •  c  In  t,  u22  =  2c(c  in  t  -  1) ,  p33  *  3c(c  In  t  -  c--^-) 
where  c  s  t/(t-l)  . 


n_ 

22 


1  as  t  -*•«*>  . 


Thus  for  large  t,  the  optimal  (o,S)  is  close  to  (0,3).  [Note:  for 

2  =  £ C - >  on  A  ,  2(a)  I  *  -1]  .  Compared  to  the  old  FGM  (1) 

*u=0 

without  iteration,  whose  max.  co-.ariance  is  attained  at  a=l  with 

2 

(max  old  Cov.)  =  (/^-Mjj)  ,  we  see  th'at 


(nax  new  Cov.) 
(max  old  Cov.) 


a  + 


.  [~>l33~u22 1 

=  [v22-yllJ 


as  t 


4.  NEGATIVE  DEPENDENCE 

LENDLS  2.  for  arbitrary,  nondegenerate  F  with  finite  mean, 
v22  ~  vil  >  y33'U22  * 

PROOF .  This  is  an  immediate  consequence  of  Theorem  1  (iir)  of 
Kadane  (1974).  According  to  the  Theorem,  there  exists  a  probability 
measure  on  the  open  interval  (0,1)  whose  expectation  is  (y^-u^)/ 

(w22"yll)- 
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A  surprising  consequence  of  the  Lemma  is:  For  the  one  iteration 
FGM  (2),  Cov(X,Y)  is  minimized  over  tl  at  (a,B)  *=  (-1,0)  regardless 
of  F  and  G  .  In  other  words,  there  exists  no  marginals  for  which 

t 

the  single  iteration  will  bring  about  a  higher  negative  covariance. 

Another  consequence  of  the  Lemma  is:  Single  iteration  can  at 
best  triple  the  positive  covariance,  as  is  the  case  in  Example  3. 

Finally,  we  conclude  with  the  observation  that,  in  view  of  (13) 

and  (25),  maximum  p  over  II  is  never  attained  at  a*B«l  as  was 

* 

tacitly  assumed  in  the  literature.  There  always  exists  (a,0)  e  A  , 
different  from  (1,1),  having  a  higher  p  value. 
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